Report for Model 2

— for both jointly and separately

2024-10-23

1. Joint Modeling

1.1 Trace plots for convergence check

The current MCMC setting is:

  • 4000000 iteration;
  • 3000000 burn-in;
  • 1000 thinning.

1.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum      1.002       1.01
## LDevsum      1.008       1.04
## dh0          1.008       1.04
## dh1          0.999       1.00
## dh2          1.011       1.05
## dh3          1.004       1.01
## dl0          1.003       1.02
## dl1          1.001       1.01
## dl2          1.009       1.04
## dl3          1.000       1.00
## 
## Multivariate psrf
## 
## 1.02

1.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

1.4 WAIC results

LevelH LevelL
DIC 1208.52394 22169.09
DIC3 1150.91711 22191.43
PWAIC 41.97006 244.77
WAIC 1177.48763 22213.35

2. Separate Modeling of High-Level

2.1 Trace plots for convergence check

The current MCMC setting is:

  • 4000000 iteration;
  • 3000000 burn-in;
  • 1000 thinning.

2.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum       1.01       1.02
## dh0           1.02       1.02
## dh1           1.01       1.03
## dh2           1.02       1.04
## dh3           1.02       1.02
## 
## Multivariate psrf
## 
## 1.01

2.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

2.4 WAIC results

H4
DIC 1347.81427
DIC3 1234.56801
PWAIC 85.41017
WAIC 1299.53354

3. Separate Modeling for Low-level

3.1 Trace plots for convergence check

The current MCMC setting is:

  • 4000000 iteration;
  • 3000000 burn-in;
  • 1000 thinning.

3.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## LDevsum       1.00       1.00
## dl0           1.00       1.01
## dl1           1.01       1.03
## dl2           1.01       1.03
## dl3           1.00       1.00
## 
## Multivariate psrf
## 
## 1.01

3.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

3.4 WAIC results

L4
DIC 22283.8370
DIC3 22301.6761
PWAIC 313.6685
WAIC 22337.1349